Wobbling kinks in 4 theory 
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We present a uniform asymptotic expansion of the wobbling kink to any order in the amplitude 
of the wobbling mode. The long-range behaviour of the radiation is described by matching the 
asymptotic expansions in the far field and near the core of the kink. The complex amplitude of the 
wobbling mode is shown to obey a simple ordinary differential equation with nonlinear damping. 
We confirm the t _1//2 -decay law for the amplitude which was previously obtained on the basis of 
energy considerations. 
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I. INTRODUCTION 

Since the early 1960s, the one-dimensional c/> 4 theory 
has been among the most heavily utilised models of sta- 
tistical mechanics and condensed-matter physics [l(. At 
the same time, it served as a testing ground for a variety 
of ideas in topological quantum field theory and cos- 
mology 3] . The equation of motion for the model reads 



(1) 



and the fundamental role in applications is played by its 
kink solution, 



4>{x,t) — tanha;. 



(2) 



The </> 4 -kinks describe domain walls in ferromagnets Q 
and ferroelectrics 0, 0] and represent elementary exci- 
tations in the corresponding partition function 7| . They 
were employed to model proton transport in hydrogen- 
bonded chains Q and charge-density waves in polymers 
and some metals 

HEt Topolo gical defects described 
by kinks are generated in large numbers during second- 
order phase transitions (Tlj : such transitions occur in a 
variety of condensed matter systems and are believed to 
have been made by different parts of the early Universe 
Q. In quantum field theory, kinks represent nonper- 
turbative classical solutions which undergo subsequent 
quantisation 113 ; one example concerns "bags" contain- 
ing fermions 1131. (For more recent quantum physics ap- 
plications see [14|.) 

Mathematically, the (f) A kink has a lot in common with 
its sine-Gordon counterpart; the two kinks are the sim- 
plest examples of topological solitons in one dimension. 
There is an important difference though; the sine-Gordon 
equation is integrable whereas the <fr theory is not. Ac- 
cordingly, the kink-antikink interaction becomes a non- 
trivial matter in the <f> 4 case @, [ID, [l||. Another (not 
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unrelated) difference is that unlike the kink of the sine- 
Gordon equation, the 4> A kink has an internal mode - 
an extra degree of freedom which allows for oscillations 
in the width of the kink. Although these oscillations 
are accompanied by the emission of radiation (another 
manifestation of the nonintegrability of the <f> A model), 
the radiation is quite weak and the oscillations are sus- 
tained over long periods of time. Since the amplitude of 
the oscillations can be fairly large, this periodically ex- 
panding and contracting kink (termed wobbling kink in 
literature, or simply wobbler) can be regarded as one of 
the fundamental nonlinear excitations of the cf> 4 theory, 
on a par with the nonoscillatory kinks and breathers. For 
small oscillation amplitudes and on short time intervals, 
the wobbler can be characterised simply as a linear per- 
turbation of the stationary kink (J2J) . However in order to 
determine the lifetime of this particle-like structure (even 
when its amplitude is small) , or characterise it when it is 
a large-amplitude excitation, one needs a self-consistent 
fully nonlinear description. 

The wobbling kink was discovered in the early numer- 
ical experiments of Getmanov (l6j who interpreted it as 
a bound state of three fundamental (i.e. nonoscillatory) 
kinks. (For a more recent series of numerical simulations, 
see [13 ■) Rice and Mele have reobtained this nonlinear 
excitation within a variational approach employing the 
width of the kink as a dynamical variable [lOL Il8| ] . Se- 
gur then constructed the quiescent (i.e. nonpropagating) 
wobbler as a regular perturbation expansion in powers of 
the oscillation amplitude . He calculated the first two 
orders of the perturbation series and noted the likely oc- 
currence of unbounded terms at the third, e 3 -, order, im- 
plying the consequent breakdown of the expansion. His 
construction was extended in Ref . [22] where the effect of 
the wobbling on the stationary component of the kink 
was evaluated. It is also appropriate to mention Ref. [23| 
where its author derived an expression for the radiation 
wave emitted by an initially nonradiating wobbler, and 
a series of publications [2J] where the interaction of the 
.agojgbler with radiation waves was studied in more detail 
and from a variety of perspectives. From the fact that 
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the energy of the wobbling kink is quadratic in the ampli- 
tude of the wobble while the second-harmonic radiation 
flux is quartic, it is straightforward to conclude that the 
amplitude decays as t' 1 / 2 [H HJ HI . 

Moving on to singular perturbation expansions, Kise- 
lev [26| studied the perturbed tfi 4 kink using the Krylov- 
Bogoliubov-Mitropolskii method. (Later, he extended 
his analysis to the </> 4 equation with a conservative au- 
tonomous perturbation [27].) His two papers have math- 
ematical rigour and a wealth of useful formulas; however 
a self-consistent system of e quat ions for the kink's pa- 
rameters was not derived in [26l . \2l\ and the long-term 
evolution of the wobbling kink has therefore remained un- 
explored. Manton and Merabet [2lJ used the Lindstedt- 
Poincare method [25[, where the expansion of the field 
is supplemented by an expansion of the frequency of the 
wobbling. Manton-Merabet's theory was successful in 
reproducing the decay law of the wobbling amplitude 
(which was previously obtained from the energy consid- 
erations [13, HH, HH). However, the Lindstedt-Poincare 
method, although efficient in finding periodic orbits, may 
lead to erroneos conclusions about nonperiodic regimes 
[25l ]. (One manifestation of this inadequacy in the case 
at hand is that the nonlinear corrections to the frequency 
become complex and time-dependent 21]; a less obvious 
difficulty is the infinite speed of the signal propagation, 
see below.) This motivates the search for new approaches 
which would be mathematically self-consistent (like the 
one in [2(| [27| ) on the one hand, and preserve the phys- 
ical insights of phenomenological expansions [U [HI on 
the other. 

The aim of the present paper is to develop a singular 
perturbation expansion of this kind. Our approach recog- 
nises the existence of a hierarchy of space and time scales 
associated with the kink+radiation system and gener- 
ates a perturbative expansion which remains uniform to 
all orders. The consistent treatment of radition requires 
also the introduction of an independent expansion of the 
far field which is then matched to the expansion near 
the core of the kink. This produces physically consistent 
and asymptotically accurate results at all space and time 
scales. In particular we will obtain a nonlinear ordinary 
differential equation obeyed by the amplitude of the wob- 
bling mode. In the follow-up paper [34] our multiscale 
approach will be used for the analysis of the wobbling 
kink driven by a resonant force. 

The basics of our method are outlined in the next sec- 
tion. In sections IIIII and IIVI we evaluate the first- and 
second-order corrections to the shape of the wobbling 
kink, and in section PVTl derive an equation for the ampli- 
tude of the wobbling mode. The asymptotic matching of 
the radiation on the short and long scale is carried out 
in sectionfVl here we show, in particular, how to account 
for finite propagation speed of radiation in a mathemat- 
ically consistent way. Finally, conclusions of this study 
are summarised in section IVIII 



II. THE METHOD 

We consider the kink moving with the velocity v. Mak- 
ing the change of variables (x, t) — > (£, r), where 

£ = x - ( v(t')dt', r = t, 
Jo 

we transform Eq.JT]) to the co- moving frame: 
v 1 — v 2 

\<t>rr - u0£ T ~ y 0$ 2 — ^ - <j> + 4> 3 = 0. (3) 

Like the authors of [28[, we shall determine the kink's 
velocity v(t) by imposing the condition that the kink be 
always centred at £ = [i.e. at x = J Q v(t')dt']. 

At first glance, the inclusion of the function v(t) is un- 
necessary: having constructed a quiescent wobbling kink, 
we could make it move at any speed simply by a Lorentz 
boost. The reason we have introduced the velocity ex- 
plicitly in Eq.([3]), is twofold. Firstly, this will allow us to 
check whether the wobbling kink can drift with a non- 
constant velocity. The soliton moving with a variable 
v(t) could obviously not be Lorentz-transformed to the 
rest frame. Secondly, we include the velocity in prepara- 
tion for the analysis of the damped-driven c/> 4 equation in 
the second part of this project [Hf. Since the damping 
and driving terms violate relativistic invariance, the ex- 
plicit introduction of the velocity becomes essential even 
when considering the damped-driven wobblers moving at 
a constant speed. 

We expand the field about the kink 4>q = tanh£: 

<p = <f>o + £01 + e 2 02 + ■ ■ ■ • (4) 

Here e is a (formal) small parameter; it will drop out of 
the final expression for the solution [see Eq. p5|) below]. 
Substituting ([4} in ([3]) and setting to zero coefficients of 
like powers of e would constitute Segur's approach which 
is expected to produce secular terms in the expansion 
[l9j . To avoid these, we introduce a sequence of stretched 
space and time variables 

X n = 6 n £, T n = e n T, n = 0,l,2,..., (5) 

which describe slower times and longer distances. In the 
limit e — ► 0, the different scales become uncoupled and 
may be treated as independent variables. We expand the 
£- and ^derivatives in terms of the scaled variables by 
using the chain rule, 

= d + ed t + e 2 d 2 + . . . , 

= D + eDi + e 2 D 2 + . . . , (6) 

OT 

where we have used the standard short-hand notation 
d n d 

° n -dX n ' Un ~dT n - 
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Substituting these expansions into the 4 equation l[3]). 
along with the series Q, and equating coefficients of like 
powers of e, we obtain a hierarchy of equations. We 
assume that the velocity of the kink is slowly varying 
and, for simplicity, that it is small, i.e. v — eV where 
V = V(T U T 2 ,...) is of order 1. 



III. LINEAR PERTURBATIONS 

At O (e 1 ), we obtain the linearisation of Eq.JT]) about 
the kink 0o = tanhAo: 



C<t> x = 0, (7) 
where we have introduced the Schrodinger operator 



C 



l£)2 



d{ - 1 + 30 2 



1«2 
2 



c^ + 2-3sech 2 A . (8) 



The general solution of the variable-coefficient Klein- 
Gordon equation ([7]) can be written as 

4>i = Cy T {X )+Ae M y w {X )+c.c. + 4> R {X ,T ), (9) 

with 

«nn 

K{p)e iu ^ To +U*(-p)e- i ^ To ] y P (X )dp. 



(10) 

Here yx and y w are eigenfunctions of the operator £ as- 
sociated with eigenvalues and I, respectively: 

y T (X ) =sech 2 A ; (11) 
y w (X ) = sech A tanh A . (12) 

The eigenfunction y w gives the spatial profile of the so- 
called internal mode, also known as the wobbling mode in 
the current context. The internal mode frequency Uq — 
y/2>. The functions y p (Xo) are solutions pertaining to the 
continuous spectrum of C: 

Cy p = ^2 + y j y p , -oo < p < oo; (13) 

these were constructed by Segur (l9j : 

3(1 - ip) 



y P (x Q ) 



JpX 



1+p 2 

3(2 -ip) 



tanh Aq ( 1 + tanh Aq ) 



(1 + tanhA ) 2 



(14a) 



4 + p 2 

We also mention an equivalent form for these solutions 



y P (x ) 



1 



JpX Q 



(1 + ip){2 + ip) 
x (2 - p 2 - 3iptanhA - 3 sech 2 A ). (14b) 



The continuous spectrum solutions are usually referred 
to as phonon modes; the phonon frequencies uj(p) are 



given by u>(p) — \J\ + p 2 > 0. Finally, the coefficients 
lZ(p) and A are complex; Co is real, and c.c. in © stands 
for the complex conjugate of the immediately preceding 
term. 

The yx eigenfunction is the translation mode; since the 
kink is assumed to be centered at Ao = 0, we let C = 0. 
Next, we will consider only localised perturbations of the 
kink, i.e. we assume that <pi — > as |Ao| — > oo. This 
means that the Fourier coefficient TZ(p) in the integral 
pop can be regarded as an ordinary function, i.e. it does 
not include any 5-function contributions. Sending To — > 
oo for the fixed finite Ao, Kelvin's formula of the method 
of stationary phase gives 

MXo,T )-> (^) ' n(0)y Q (X )e 2 * T °+^ 4 + c.c. 

Therefore 4>r is a slowly-decaying wavepacket which will 
be dominated by the zero-wavenumber radiation after 
other harmonics have dispersed away. 

Since we are interested in the evolution of the wobbling 
mode and not that of a general localised initial condition, 
we set TZ(p) = 0. Therefore the first-order perturbation 
is taken in the form 

0i = A(Xi, . . . ; Ti, . . .) sech A tanh X e iuJoT ° + c.c. 

The amplitude of the wobbling mode, A, is constant with 
respect to Ao and To but may depend on slower times 
and longer distances. 



IV. QUADRATIC CORRECTIONS 

At the second order in the perturbation expansion we 
arrive at a nonhomogeneous variable-coefficient Klein- 
Gordon equation 

^Dlfc+Cfo =F 2 (A ,...;T ,...), (16) 
where the forcing term is 

F 2 = (dodi - A>£>i)0i " 30o0 2 + VD d (f>i 

+^£>il/d o 0o ~ lv 2 d%<p - (17a) 
Substituting for <po and <pi, this becomes 



F 2 = -6\A\ 2 sech 2 A tanh 3 A + -D X V sech 2 A 



+V 2 sech 2 A tanh A + [di A{2 sech 3 A 
— sech Ao) — iuioDiA sech Ao tanh Ao 
-iuj VA(2 sech 3 A - sech A )] e 



+ C.C. 



-3A 2 sech 2 A tanh 3 Aoe 2 ^ + c.c. (17b) 



The T -independent term in Eq. (|17b[) and the term pro- 
portional to e iuJ ° T ° are resonant with the two discrete 
eigenmodes of the operator in the left-hand side of (|16p , 
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while the term proportional to e 2luJ ° To is resonant with 
its continuous spectrum. The latter part of the forcing 
is localised in the region near the origin and acts as a 
source of radiation which spreads outward from there. 

We discard the homogeneous solution of Eq. (Tlo]) for 
the same reason as we have discarded most terms in the 
solution of Eq. {7J ; namely, we do not want the evolution 
of the wobbling mode to be shaded by dispersive tran- 
sients. Hence the solution that is of interest to us will 
consist only of the harmonics present in Eq. pT)) : 

- ^ 2 0) + ip^e M + c.c. + <p^e 2 ^> T <> + c.c, (18) 

where (p 2 °\ <P;P an d are functions of Xq which satisfy 
the three linear nonhomogeneous equations: 

Clp [ ° ] = -Q\A\ 2 scch 2 X tanh 3 X Q 
+ iL>il/sech 2 X Q + V 2 sech 2 X tanhX , (19) 

(£ - |)^ X) = <9iv4(2sech 3 X - sech A ) 
— ioJoDiA sech X tanh X 
+iuj VA(2 sech 3 X„ - sech A ), (20) 

and 

(C - 6)^4 2) = -3A 2 sech 2 X tanh 3 X . (21) 

(The functions (f 2 °\ tf-P anc ^ Pp can a ^ so depend, para- 
metrically, on Xi,X 2 , ... and T l7 T 2 , ....) 

The homogeneous solutions of the first two of these 
equations are given by the eigenfunctions of the operator 
£, Eqs. fTTj) and (JT3J) . According to the Fredholm alter- 
native, the nonhomogeneous equations admit bounded 
solutions if and only if their right-hand sides are orthog- 
onal to the corresponding homogeneous solutions. For 
this to be the case, we must set D{V = and D\A = 0. 
The variation of parameters yields then 

ipf ] =2\A\ 2 sech 2 X Q tanh A + 

(V 2 \ (22) 

+ f — -3|A| 2 J A sech 2 A 

and 

ip^ = -(diA + iu VA)XosechX tanhX . (23) 



whence 

A = A(X 2 ,X 3 ,...;T 2 ,T 3 ,...)e-^ oVX K (25) 

We also note the terms proportional to Xq sech 2 Xq 
in Eq. ([22"l) . These terms do not grow bigger than the 
previous term, 4>q — tanhXp, yet they become larger 
than the difference 4>q — 1 as Xq — > oo and <j)Q + 1 as 
Xq —t—oo. If we attempted to construct the asymptotic 
expansion of the function <p— 1 at the right infinity or the 
function <f> + 1 at the left infinity, the terms in question 
would cause nonuniformity of these expansions. Since 
the function Xq sech 2 Xq is nothing but the derivative 
of tanh(fcAo) with respect to fc, these terms represent 
the variation of the kink's width. Hence the potential 
nonuniformity of the expansion can be avoided simply 
by incorporating them in the variable width [see Eq. (|4"5|) 
below] . 

We now turn to the remaining nonhomogeneous equa- 
tion, Eq. (|2ip . The variation of parameters gives 

v4 2) =A 2 A(Xo), (26) 

with 

fi (X ) = | { 6 tanh X sech 2 X Q 

+ (2 + ik tanhXo + sech 2 X )[J 2 *(Jf ) - J™]e %k ° x ° 

+ (2- i/c tanhX +sech 2 Xo)J 2 {Xo)e~ lkoXo }. (27) 

Here the function J 2 (X ) is defined by the integral 

e lkoi sech 2 f d£ (28) 

-oo 

with k = i/8. The constant J%° is the asymptotic value 
of J 2 (X ) as Xq — > oo: 

J 2 °°= lim J 2 {Xq). (29) 

Xq — >CO 

The two constants of integration were chosen such that 
the solution (|26 [) - (f5T|) describes right-moving radiation 
for positive Xq and left-moving radiation for negative 
Xq . It is not difficult to show that f\ is an odd function; 
we will use this fact in what follows. 



V. RADIATION IN THE FAR FIELD 



Although the function ip 2 ' decays to zero as |X | — > 
oo, the product eip^ becomes greater than the first-order 
perturbation yi(Xo) for each fixed e and sufficiently large 
|Ao|. Consequently, the term e 2 (f> 2 in the expansion (j4|) 
becomes greater than the previous term, e <f>i, leading 
to nonuniformity of the expansion. In order to obtain 
a uniform expansion, we set this "quasisecular" term to 
zero: 

dxA + iuj VA = 0, (24) 



The function (|2"6")l is bounded but does not decay to 
zero as \Xq\ — > oo. This fact presents a problem, both 
for the consistency of our method and for the physical 
interpretation of the resulting solution. Mathematically, 
the term e 2 (f> 2 turns out to be greater than the previ- 
ous term in the expansion ((4]) for sufficiently large \Xq\. 

As we have mentioned in connection with the term ip 2 , 
this leads to nonuniformity of the expansion. Physically, 
the problem is that any variation of the amplitude of the 
wobbling mode, A, on the time scale T 2 , will result in 



5 



a simultaneous change in the amplitude of the radiation 
tail for all values of Xq, from the origin to the plus- and 
minus-infinity. This is obviously in contradiction with 
the finiteness of the velocity of signal propagation in a 
relativistic theory [which is bounded by 1 in the dimen- 
sionless units of Eq.JT]).] 

The problem stems from the fact that the equation (TIB)) 
and, therefore, equation (p?Tj) , were obtained under the as- 
sumption that, in the expansion (j4]), the second term is 
smaller than the first one, the third one is smaller than 
the second, and so on — more precisely, that e<px/4>o — > 0, 
£ 2 02/(e0i) — * 0, and so on, as e — * 0. This assumption 
turns out to be only valid on the short scale and there- 
fore, the equation (f2Tj) is only meant to hold for distances 
X = 0(1) but not X = Oie- 1 ) or longer. The interval 
of Xq where e0„+i/0 n — ► as e — ► will be referred to 
as the "inner" region in what follows. Eqs. (fT6f and Q21p 
are therefore valid in the inner region. 

To obtain a uniform expansion on the whole axis, we 
also consider two "outer" regions — one with X > and 
the other one with X < 0. We define the outer regions 
by the requirement that |X | be greater than ^lne -1 . 
Note that the outer regions overlap with the inner region. 
For example, the values Xq = ±| In e _1 are clearly in the 
outer regions; on the other hand, we have e0i/0o — * 0, 
e 2 02 1 - ► 0, etc. for these Xq and so they belong to 
the inner region as well. 

In the right outer region, we expand in the power 
series 

0=1 + e 2 (j)2 + e 4 4 + ■ • ■ , (30a) 

and in the left outer region, we let 

= -l + e 2 2 + e 4 4 + .... (30b) 

Substituting these, together with the expansions ([6]), in 
Eq.©, the order e 2 gives 

±£ 2 2 + £0 2 = 0, 

where £ = — \8q + 2 is the far-field asymptotic form of 
the operator ©. The solutions of this equation in the 
right and left outer regions are, respectively, 

2 = j B+e l ^+ T °- k + x ^ + c.c. (31a) 

and 

02 = -JB-e i{u} - T °- k - Xa) + c.c, (31b) 

where lj 2 . = k\ +4, and the amplitudes B± are functions 
of the "slow" variables: B± = B±(X t , T x , ...). The 
normalisation constant J will be chosen at a later stage, 
and the negative sign in front of B- is also introduced 
for later convenience. 

Eqs. ([3i"l) should be matched to the solution in the inner 
region, Eg. 1(15)1 with coefficients as in $2$ , <$5§ , and ((2"rJj) . 
To this end, we take the values X — ±| In e _1 (which, as 
we remember, belong to the overlap regions). For these 



X Q , we have \X X \ ^0(e\ne- 1 ), \X 2 \ = O (e 2 In e^ 1 ) 
and so X\ — > 0, X 2 — * 0,..., as e — > 0. The solutions (pIT)) 
become, in this limit: 

02 = ±JB ± (0,0, ...;T 1 ,T 2 , ...)e l{uJ±To ~ k±Xa) + c.c. 

On the other hand, letting \X \ = line" 1 and sending 
e -> in Eqs.©, and |(2BJ|, we get 

02 =±(2-^ )J 2 °°A 2 (0,0,...;T 2 ,r 3 ,...)e i ( 2 ^ To ^oXo) +ca; 

where the top and bottom sign pertain to the positive and 
negative Xq, respectively. Choosing J = (2—iko)J%° and 
equating the above two expressions, we obtain w± — 2u> , 
k± = ±feo, and 

B±(0,0,...;Ti,T 2 ,...) = ^ 2 (0,0,...;T 2 ,T 3 ,...). (32) 

Eqs. ((52^) can be regarded as the boundary conditions 
for the amplitude fields B + and B_ . Equations govern- 
ing the evolution of these functions of slow variables can 
be derived at higher orders of the (outer) perturbation 
expansion. Namely, the solvability condition at the order 
e 3 yields 

(dodi - D D 1 + VdoDa)^ = 0. (33) 
Substituting from (pH) . this becomes 

DiB± + ^diB± + ik ± VB± = 0, (34) 

whence 

B ± = e- 2iu ° vx *B±(X 1 ,X 2 ,...;T 1 ,T 2 ,...), (35) 

where B± satisfy a pair of linear transport equations 

DiB+ + c Q d 1 B+ = 0, Xi > 0, (36a) 
D X B- - codxB- =0, Xx < 0, (36b) 

with co = ko/(2u>o)' Note that cq is nothing but the 
group velocity of the radiation waves with the wavenum- 
ber k,Q\ cq = (dui/dk)\ k _ k , where to = V4 + k 2 . 

Solution of equations (f36| with the boundary condition 
((52")) is a textbook exercise. Assume that the functions 
B± satisfy the initial conditions 8^(^,0) = B^iX^ 
(for X x > 0) and B-(X u 0) = B^\X X ) (for X x < 0), 
with some function B^°\Xi) defined on the whole axis 
-oo < X x < oo, with B^(Xx) -»• as |Xi| -> oo. 
(We have suppressed the dependence on the variables 
X 2 , X3, T2, T3, ... for notational convenience.) In the 
region Xx > coTx, the solution to the equation (|36a[) 
with the above initial condition is given by B+(Xx, Tx) — 
B^ (Xx—cqTx). This solution represents an envelope of a 
group of second-harmonic radiation waves, moving to the 
right with the velocity cq. Importantly, the amplitude B+ 
in this region is not related to the wobbling amplitude A 
and so no information from the core of the kink can reach 
this region. In the region < Xx < cqTx, the solution 
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to Eq. (|36aP is determined by the boundary condition in- 
stead: B + (X 1 ,T 1 ) = A 2 (0;0). This result implies that 
the moving envelope has the form of a propagating front, 
leaving Z5+ flat and stationary in its wake. In a similar 
way, on the negative semiaxis we have a front moving 
with the velocity — Co and leaving 0_(Xi,Ti) equal to 
the constant A 2 (0;0) in its wake. 

The above analysis has two shortcomings. One draw- 
back is that we have restricted ourselves to groups of radi- 
ation waves with the characteristic length and time scale 
of order e^ 1 . A natural question therefore is whether 
variations with larger space and time scales (e.g. vari- 
ations on X 2 and T 2 scales) could not propagate faster 
than Co- Another latent defect is that the solutions for 
B±(Xi,Ti) that we have constructed, will generally be 
discontinuous along the lines X\ = ±c Ti. To address 
both of these issues, we proceeed to the order e 4 of the 
outer expansion where the solvability condition for the 
second harmonic gives 



i(2lu D 2 + k 8 2 )B ± + \{D\- d 2 )B± 
-iVfaDx ~ 2^ 9i)B ± - \v 2 k 2 B ± = 0. 



Eliminating D\B± using this becomes 



VI. DECAY LAW FOR THE WOBBLING 
AMPLITUDE 

Returning to the original, "inner" , expansion ^ and 
collecting terms of order e 3 gives the equation 



\E%fa + £0 3 = ^3 



(39a) 



where 



F 3 = (do&i - A)^i)02 + {d a d 2 - D D 2 )fa 



m-Di) 



4>\ - 600^102 + VD d Q (j> 2 



+ VDM! + md o 0i + - |y 2 ^0i. 

(39b) 

Having evaluated F3 using the known functions 0o, 0i 
and 02, we decompose the solution 03 into simple har- 
monics as we did at O (e 2 ) . The solvability condition for 
the zeroth harmonic in equation (|39b[) gives D 2 V = 0, 
which means that V remains constant up to times t ~ 
e~ 3 . The solvability condition for the first harmonic pro- 



duces 



where 



i^D 2 A + t\A\ 2 A-V 2 A = 1 (40) 



iD 2 B± ± ic d 2 B ± - iVd x B± - ^±d\B ± = 0, (37) 



where u kk = (d 2 u> / dk 2 )\ ko = (4w 2 — k 2 )/(8ujQ) is the 
dispersion of the group velocity of the radiation waves. 
Combining Eq. (|37f with (|34|) . we obtain a pair of equa- 
tions in the original space and time variables: 



id t B± ± ic Q d x B± t vk B± - ^d 2 B ± = 0. (38) 



The pair of linear Schrodinger equations (|38|) govern 
the evolution of the radiation amplitudes over times and 
distances as large as e -2 ; if we want to have a description 
on even a larger scale, we simply need to include equa- 
tions from higher orders of the outer expansion. Solutions 
of Eqs. ([55]) with the boundary conditions B± = A 2 at 
x = vt and B± = at x = ±00 have the form of slowly 
dispersing fronts propagating at the velocities ±co and 
interpolating, continuously, between A 2 and 0. As in our 
previous description exploiting the transport equations 
(I36[) and valid on a shorter space-time scale, perturba- 
tions of A 2 cannot travel faster than cq, the group veloc- 
ity of radiation. 

Thus, by introducing the long-range variables B±, "un- 
tied" from the short-range amplitude A, we have restored 
the finiteness of the velocity of the radiation wave prop- 
agation. By introducing the outer expansions, we have 
also prevented the breakdown of the asymptotic expan- 
sion at large distances. 



C = 6 / sech 2 X tanh 3 X 



■ sech Xq tanh Xq 



-3X sech 2 X Q + A(A )j dX . (41) 

Out of the real and imaginary part of £, the imaginary 
part is more important; it can be easily evaluated ana- 
lytically: 







37T 2 fc 



sinli (nk /2) 



= 0.04636. 



The real part was computed numerically: 
( R = -0.8509. 



(42) 



(43) 



Denoting eA = a the "natural" (unsealed) amplitude 
of the wobbling mode, and recalling that v — eV and 
A t — e 2 D 2 A + (e 3 ), we express the amplitude equation 
(]40[) in terms of the original variables: 



ia t = — \a\ a + —v a + O [\a\ 



(44) 



Eq. (j44j) contains solvability conditions at all orders cov- 
ered so far — they arise simply by expanding the deriva- 
tive d/dt as in Eq.©. Unlike the amplitude equation 
D\A — which only governs the evolution for times 
t ~ e _1 , and unlike the equation fll))) which only holds 
on the timescale t ~ e -2 , the "master equation" (|4"4")l is 
applicable for all times, from t — to t ~ eT 2 . 
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FIG. 1: The decay of the free wobbling kink. Crosses: 
\a(t)\ as obtained from the direct numerical simulations of 
Eq.Q with the initial conditions in the form cj> = tanh 2 + 
2ao sech x tanh x, (j>t = 0, with ao = 0.05. Solid curve: equa- 
tion |g6| with | a (0)| = 0.0492. The inset shows the first 250 
periods of oscillation; during this time the amplitude drops 
by less than 10%. 



tions to the frequency of the wobbling: 

u=u J0 [l-y + ±( R \a\ 2 + O(\a\ i )}, (47) 

with < as in (j43|) . (Note that ui is the frequency of 
oscillation of the "full" field cf>, not just of the amplitude 
a.) The |a| 2 -term here is a nonlinear frequency shift from 
the linear frequency u>o — V3; as time advances, this 
term decays, slowly, to zero. The v 2 -term comes from 
the transverse Doppler effect. We could have obtained 
this term simply by calculating the wobbling frequency 
in the rest frame and then multiplying the result by the 
relativistic time-dilation factor \/l — v 2 (which becomes 
1 — \v 2 for small v). 



VII. CONCLUDING REMARKS 



The master equation (|44[) is the final result of the 
asymptotic analysis. All the conclusions about the be- 
haviour of the wobbler's amplitude shall be made on the 
basis of this equation. We could extend the range of 
applicability of the master equation beyond times of or- 
der e -2 by continuing our perturbation analysis to higher 
orders of e. However, corrections to the equation (f4"4")) ob- 
tained in this way would be smaller than the terms that 
are already in the right-hand side of (|4~4"1) and would not 
affect our conclusions based on (|44|) in its present form. 

The absolute value of a is governed by the equation 

||a| 2 =^ C/|a| 4 + O(| a | 6 ). (45) 

Previously this eq uation was obtained using heuristic 
considerations [2Ct l2ll . |23| . Since Ci > 0, the ampli- 
tude of the wobbling is monotonically decreasing with 
time: a constant emission of radiation damps the wob- 
bler. Dropping the O (|a| 6 ) correction term from (|45| . 
the decay law is straightforward: 

, m|2 = KQ)I 2 = KQ)I 2 

1 Ul 1 + lu (i H0)\H 1 -I- 0.08030 x |a(0)| 2 i' 

(46) 

When a(0) is small, the decay becomes appreciable only 
after long times t ~ |a(0)|~ 2 . The decay is slow; for times 
t > 12.5 x |a(0)|- 2 , Eq.fjll gives |a| ~ r 1 / 2 . 

We have verified the above decay law in direct numer- 
ical simulations of the full partial differential equation 
([T]). (The details of our numerical algorithm have been 
relegated to the Appendix.) As the initial conditions, we 
took cj>{x, 0) = tanhx + 2ao sech x tanh x with some real 
ao and (f>t(x,0) — 0. After a short initial transient, the 
solution was seen to settle to the curve (|46| with |a(0)| 
close to ao, see FigfTJ 

The equation l|44p gives us the leading-order contribu- 



In this paper, we have formulated a singular perturba- 
tion expansion for the wobbling kink of the </> 4 model. 
Unlike the previously published singular perturbation 
theories based on the Krylov-Bogoliubov and Lindstedt 
methods, our approach exploits the existence of multi- 
ple space and time scales in the kink+radiation system. 
Some aspects of our scheme are standard to the method 
of multiple scales; some other ones (e.g. the appearance 
of the quasisecular terms) are less traditional. We partic- 
ularly emphasise our novel treatment of the long-range 
radiation and the infinite propagation speed paradox. 
The final result of the asymptotic analysis is the am- 
plitude equation for the wobbling mode, Eq. p^)) . Using 
this equation, we evaluate the nonlinear frequency shift 
and decay rate of the wobbler. 

The coupling of a spatially localised temporally pe- 
riodic excitation to radiation modes via a nonlinearity 
was discussed previously in several contexts. In particu- 
lar, Ref. described the decay of the internal mode of 
the nonlinear Schrodinger soliton, in the equation with 
a general nonlinearity. (For rigorous estimates, see e.g. 
(30j.) In Ref.[3l|, the dynamics of the soliton's internal 
mode was considered in the nonlinear Schrodinger equa- 
tion with the parametric forcing and damping. Next, the 
authors of Ref . [32] studied the persistence of a localised 
linear impurity mode in the cubic Klein-Gordon equa- 
tion. [We note that although our Eq.([T]) can also be cast 
in the form of an equation with an impurity potential 
- by letting <f) — + X — the resulting Klein-Gordon 
equation satisfied by x does not fall into the class of sys- 
tems covered by the analysis in that paper.] We also 
mention an earlier article [33| where a similar problem 
was considered for the nonlinear wave equation. 

We conclude our study by producing the perturbation 
expansion of the wobbling kink in terms of the original 
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variables: 



4>{x, t) = tanh 



l-3|o| 



+ asech^tanh^ n(t " 1 ' ?) + c.c. 

+2|a| 2 sech 2 etanhe+a 2 /i(0e 2M0(t " t,?) +c.c.+O(|a| 3 ) . 

(48) 

Here £ = x — vt; the complex function a(t) satisfies an 
ordinary differential equation (j44[) , and f± (£) is given by 
Eq. ({2"T)) . Note that we have incorporated two Xo sech 2 Xq 
terms of the sum (|2"2"|) into the variable width of the 
kink. The expansion ([48]) is only valid at the length scale 
|£| = O (1); for larger distances one has to use the outer 
expansions (|30p with coefficients determined in section 

The first term in (|48)l describes a moving nonoscillatory 



kink with the width decreasing (to the value of yl — v 2 ) 
on the timescale t ~ | a, | 2 . The second term describes 
the wobbling mode; the third gives the quasistationary 
correction to the shape of the kink induced by the wob- 
bling, and the last term accounts for the second-harmonic 
radiation from the wobbler. 

The first term in (|48[) is manifestly Lorentz-covariant. 
The other terms can also be cast in the relativistically- 
covariant form if we replace £ with £/vl — v 2 in sech£ 
and tanh£ (this is correct to the order of v 2 ), and write 
ae 



iu>0 (*—«£) 



a exp 



iw (l + iCfll«| 2 ) 



Here we used Eq. (|4T|) and neglected terms of order |a| 4 x 
and \a\ 5 t. (We remind the reader that v and |o| are con- 
sidered to be small quantities, of the same order of small- 
ness.) 
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APPENDIX A: NUMERICAL METHOD 

In this Appendix we outline some relevant details of 
the numerical method that we used to verify predictions 
of our theoretical analysis. 

Equation JT]) was simulated using an explicit finite- 
difference scheme on a grid of step size Ax = 0.1 and 
At = 0.05. The simulations were performed on the in- 
terval —L<x<L, where L was chosen large enough to 
prevent the kink exiting the domain of integration. (Typ- 
ical values of L were of the order of 1000.) We imposed 
the free-end boundary conditions. 

In order to prevent the radiation reflecting back from 
the boundaries of the system, damping was introduced 
near the edges to absorb the radiation. That is, we added 
to the </> 4 equation an absorbing term 7(2;)^, with 



7(2) = 



x-(L-WQ) 
100 

z+(L-100) 



100 



for x > L - 100; 

for x < -L + 100; 
otherwise. 



The position xo(i) of the wobbling kink was deter- 
mined from the location of the zero crossing. The am- 
plitude of the wobbling mode was measured by tak- 
ing the profile cf>(x,t), subtracting the reference kink 
tanh[ir — xo(t)], and assuming the odd component of 
what remains to be the first-harmonic wobbling mode, 
asechAotanhAoe 1 " 01 " + c.c.. This technique, of course, 
furnishes only a first-order approximation to the ampli- 
tude because of the higher order terms in the perturba- 
tion expansion. Interpolation and smoothing were ap- 
plied to counter the effects of the discreteness of the x 
values and the various oscillations occurring on the fast 
time scale. 
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